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The use of coherent x-ray beams has been greatly developing for the past decades 1 . 
They are now used by a wide scientific community to study biological materials 2 , 
phase transitions in crystalline materials 3 , soft matter 4 , magnetism 5 " 6 , strained 
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structures " , or nano-objects . Different kinds of measurements can be carried 
out: x-ray photon correlation spectroscopy 10 allowing studying dynamics in soft 11 
and hard 12 matter, and coherent diffraction imaging 3 enabling to reconstruct the 
shape 14 and strain 15 of some objects by using methods such as holography 16 or 
ptychography 17 . In this article, we show that coherent x-ray diffraction 18 (CXRD) 
brings a new insight in another scientific field: the detection of single phase defects 
in bulk materials. Extended phase objects such as dislocations embedded in the 
bulk are usually probed by electron microscopy or X-ray topography. However, 
electron microscopy is restricted to thin samples, and x-ray topography is 
resolution-limited. We show here that CXRD brings much more accurate 
information about dislocation lines (DLs) in bulk samples and opens a route for a 
better understanding of the fine structure of the core of bulk dislocations. 

Understanding the origin and behaviour of topological defects in hard condensed 
matter is a challenging problem. Their study has always been of prime importance 
because of the large influence they have on the physical properties of materials 19 . Good 
examples are the plasticity of crystalline materials, determined by the mobility of 
dislocations, the electronic conductivity, strongly affected by the density of dislocations, 
or even more complex phenomena like the sliding motion of charge/spin density waves 
(CDW/SDW) mainly related to the dynamics of climbing DW dislocations. Only a few 
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techniques, such as x-ray topography or electron microscopy are able to probe the 
properties and dynamics of dislocations. We show that CXRD not only is a new tool to 
image DLs located deep inside samples, but also provides complementary information 
to the so far existing techniques. Another strength of CXRD is its ability to probe phase 
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defects of any x-ray diffractive order parameter. Its sensitivity to CDW and SDW 
dislocations has already been reported. In addition, combination of CXRD with 
topography proves efficient to get rich and accurate information of extended bulk 
dislocations and allows to directly image them. 

The combination of CXRD and topography takes advantage of the behaviour of a 
light beam limited in size by slits. If an aperture a defines the size of a ^-wavelength 
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monochromatic beam, the latter keeps parallel up to a distance dFF=a I AX, and then 
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diverges linearly. If both sample and detector are inserted in the parallel beam regime, 
i.e. below a distance dFF from the slits, a topography image is obtained. In practice, this 
corresponds to large apertures a. On the contrary, a CXRD image is obtained when the 
sample receives a parallel beam while the detector receives a diverging beam i.e. if dFF 
is intermediate between sample and detector distances from the slits, and when the 
aperture a is lower than the transverse coherence length of the beam. The two setups 
are displayed in Fig.l (see section Methods for further details about the setups). 




Figure 1 | Coherence and topography setups used for the experiment. A 

monochromatic beam was selected by a double-crystal Si(1 11) monochromator. 
No focusing has been used to get the required parallel beam for topography. 
Source slits were used as secondary source. For topography, attenuators were 
also introduced. For coherence, slits were used before the sample to select only 
the coherent part of the beam. See section Methods for further information. 



To demonstrate the strength of combining CXRD and topography, we present a 
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study of the well-known oxygen-induced dislocations in silicon . This sample 
processing stabilizes two kinds of defects: stacking faults (SF) bounded by partial 
dislocation loops (Burgers vectors bs=l/3<l 1 1>), and perfect dislocation loops, with 
Burgers vectors bp=l/2<110>. As shown in Fig. 2a, the latter are dissociated into two 
partial dislocations, with Burgers vectors bi and b2 respectively, so that bp = bi + b2. 
The partial dislocations are called glissile if their direction and Burgers vector belong to 
the same {111} glide plane, as sketched in Fig. 2a, or sessile, if not. As a consequence, 



2 



a SF stands between the two partial DLs provided one at least is glissile. Topography 
is first used to locate the defects to be studied. As the dislocation loops are extended 
over 100 \im on the average, the several hundred micrometers size beam used in 
topography is very convenient to get a rapid overview of the defects positions and 
concentration. A 5 jam beam is then used to make a CXRD study of the selected defect. 




Figure 2 | Coherent Bragg reflection associated to a dissociated 
dislocation, a, Dissociation of a perfect b=1/2[1 10] dislocation developing 
along a [-112] direction into two glissile partials: bi=1/6[121] and b 2 =1/6[21-1]. 
The yellow plane is a stacking fault. Superimposed to the schematic atomic 
lattice, the colour scale refers to the phase shift induced by the presence of the 
dissociated dislocation; b, measured (220) Bragg reflection in the presence of a 
dissociated dislocation. A splitting due to the phase shift is observed along with 
a streak perpendicular to the stacking fault; c, Calculated diffraction pattern 
associated to a dissociated dislocation, in good agreement with the measured 
intensity distribution shown in b. 
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Thanks to its high sensitivity to phase shifts compared to standard diffraction, 
CXRD is a very powerful technique to probe dislocations. The n phase shift induced by 
a dislocation is expected to lead to a destructive interference at the Bragg position. In 
the case of a dislocation located in the middle of the coherently illuminated volume, the 
Bragg peak is expected to split into two parts, each of them having the same angular 
width as that of a Bragg reflection associated to a perfect lattice, and approximately half 
its maximum intensity (for more details, see Supplementary Discussion 1 in 
Supplementary Information). This dislocation induced splitting has been experimentally 
recorded on the (220) Bragg reflection, as shown in Fig. 2b. In addition, despite the 
thinness of the nm-size SFs due to dissociation, the data reveal streaks appearing 
together with the splitting, along <111> directions, i.e. perpendicular to the SF. A 
calculation of the expected (220) reflection in presence of a dissociated dislocation 
shows a very good agreement with the measured intensity distribution (Fig. 2c). 
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Figure 3 | Principle of the scanning CXRD. a, Typical coherent diffraction 
pattern obtained on the (220) reflection in presence of a dissociated DL, with a 
10 mdeg detuning of the sample. The truncation rod is split (1) due to the phase 
shift induced by the dislocation, and a diffuse streak (2) appears along the 
splitting, perpendicular to the SF located in between the dissociated lines, b, 
The typical patterns shown in (a) only appear in the close vicinity of the DL. The 
diffuse streak is always perpendicular to the DL (and SF) direction. The red 
lines correspond to DLs along the [-110] direction (sessile), and the blue ones 
along [-11 ±2] directions (glissile). c, The measured DLs directions are in 
agreement with the expected directions and angles for dissociated dislocations 
(a = 109.47°, p= 125.26°). 
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This unique signature under a coherent light allows imaging such extended 
defects by performing a spatial scan across the area of the dislocation, and taking 
advantage of the difference in intensity at the Bragg position when the DL is or is not in 
the illuminated volume. When a perfect region is probed, the reflection is unique, and a 
maximum appears, whereas when the DL is in the illuminated volume, the peak is split, 
and the value at the Bragg position is lowered. An enhanced contrast is obtained by 
slightly detuning the sample to avoid the intense Bragg maximum. The measured 
distribution of intensity is then composed of a cut through the truncation rods (TRs) 
associated to the main (220) Bragg reflection, and to the diffracting slit fringes (for 
more details, see Supplementary Discussion 2 in Supplementary Information). In that 
situation, Bragg reflection splitting is observed on the proper associated TR (label 1 in 
Fig. 3 a) and the signal from the diffuse streak becomes thereby more contrasted (label 2 
in Fig. 3a). Fig. 3b displays a typical 2D scan performed with a 5 um coherent beam 
over a region beforehand selected with topography. As expected, the splitting and 
streaks only appear in the close vicinity of the DLs, which is the key to perform imaging 
of the line with beam-size resolution. The DLs directions are determined by delimiting 
in space the regions where the splitting and streaks occur, and are found to be the 
expected ones for a perfect dislocation loop with b p =l/2[l 10] (see Fig. 3c). Moreover, 
the streaks appear systematically perpendicular to the partial dislocations possible 
{111} glide planes. The same method was repeated across another dislocation loop (see 
Fig. 4a), first selected with topography (Fig. 4b). A 200 x 200 um 2 region was scanned 
making steps having the same size as the beam footprint in each direction. The blue 
valley of Bragg peak minima in Fig. 4c. directly gives the position of the dislocation 
loop lines and is in good agreement with the topography image. However, the present 
technique allows going beyond topography, as it informs about the presence of a 
dissociation of the DL if any, through the local measurement of additional streaks in the 
diffraction pattern. 




Figure 4 | Image of a dislocation loop obtained with scanning CXRD. a, 

Coherent diffraction patterns obtained in the area of a dislocation loop, 
beforehand selected making b, a topography image of the same dislocation 
loop, c, Imaging performed taking the intensity at Bragg reflection position for 
each pattern shown in a. When the DL is in the probed volume, the value is 
minimum, which allows following it in space. The perfectly ordered regions 
appear in red, the DL in blue. The topography and CXRD images are in good 
agreement. 
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In conclusion, this article proves the high sensitivity of CXRD to any dislocation 
embedded within a few micrometers in the bulk, and shows that direct imaging of bulk 
dislocations by CXRD is possible. The combination of CXRD and topography is a very 
flexible method which can be used at every beamline dedicated to diffraction, to first 
locate the object in the real space and second to probe its fine structure through its 
Fourier transform in the reciprocal space. 

In the canonical example of dislocation loops in silicon, the obtained CXRD 
patterns display a destructive interference at the Bragg position as well as a diffuse 
streak emanating from the fine structure of the dislocation core region. The diffraction 
pattern is compatible with the presence of dissociated dislocations, as observed by 
electron microscopy studies in thin samples. Only the destructive interference allows 
performing the imaging of the DL itself, which highlights the need of a coherent beam 
for this study. 

To push the study further, such improvements as 3D data collection, the use of 
sub-micrometer beams 4,25 and of reconstruction algorithms e.g. ptychography would 
give access to unprecedented resolution and to the strain field associated to the 
dislocation line. Besides, it also represents a unique opportunity to probe the core of 
topological objects embedded in the bulk, with an unequalled resolution and in a non 
destructive way. In that case, a more detailed model would have to be used to describe 
the dislocation core with more accuracy, although in the present case the simple phase 
model nicely explains the different features of the diffraction pattern. 

With the quest for new materials like SiC and the advent of new processes like 
bonding in view of improving device performances, the study of topological defects, in 
particular in semiconductors, is experiencing a revival for which this work opens 
significant perspectives. This approach also benefits from the easiness to use various 
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sample environments, such as high pressure , AFM , magnetic/electric fields or 
Laser pulses 29 . For these in-situ experiments, the new XFEL sources may enable to 
probe fast dynamics of dislocations at the femtosecond timescale. 



METHODS 
Experimental details 

Preliminary experiments were performed at the BM02 beamline of the ESRF. The data 
presented here were obtained during an experiment performed at the CRISTAL 
beamline of the synchrotron-SOLEIL. The source is a U20 undulator, and a Si(lll) 
double-crystal monochromator (DCM) selected a wavelength ^=1.7462 A, with 
AMk=\A 10" 4 . Focusing devices (DCM sagittal bender and focusing mirrors) were not 
used in order to have a parallel beam for topography. Rejection of harmonics was 
ensured by detuning the DCM at about 50% of the maximum value of the rocking 
curve. Secondary slits were used for beam collimation 13 m upstream of the sample and 
attenuators prevented saturation of the detector. The coherent beam was generated by 
setting the secondary slits to 200 um, leading to a transverse coherence length of 11.3 
um at the sample position. Detection was performed with a 13 um pixel size ANDOR 
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direct illumination CCD camera standing 2.2 m after the sample. Slits located 8 cm 
upstream of the sample were closed at 7 um, allowing to get a parallel beam at the 
sample position, and a divergent one at the detector position. The total transverse 
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coherence degree was then about 80%. 
Sample preparation 

The 18x9x0.5 mm silicon single crystal was grown by the Czochralski method, and 
then annealed during 35 hours at 1100 °C under an oxygen atmosphere. The 
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concentration in oxygen impurities was then 10 atoms/cm . The surface was oriented 
(110). A 4 mdeg mosaicity was measured by performing rocking scans on the (220) 
reflection. The sample was provided by G. Rolland from CEA-LETI Grenoble. 
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Supplementary Discussion 1 
Destructive interference due to the presence of a dislocation line in a 

coherently illuminated volume 

The problem can be first simplified to a one dimensional model. If one considers a 
perfectly ordered lattice (blue line on Supplementary Figure l.a), the associated phase is 
constant in space, and can be arbitrarily chosen to be (blue line on Supplementary 
Figure l.b). Any Bragg reflection associated to this lattice is unique, and its width is 
inversely proportional to the dimension of the system (blue curve on Supplementary 
Figure l.c). By introducing a 71 phase shift in this lattice, one obtains the lattice and its 
phase variations shown with red lines on Supplementary Figure l.a and b. The Bragg 
reflection associated to this lattice containing a single n phase defect is represented in 
red on Supplementary Figure l.c. The intensity drops down to zero exactly at the 
expected Bragg position. This is due to the destructive interference that occurs between 
the amplitudes scattered by the two out-of-phase volumes. Let us consider that the 
lattice has a length L, and that Ai(q) is the amplitude scattered by the volume located 
between x = and x = L/2, and A2(q) the amplitude scattered by the volume located 
between x = L/2 and x = L. The total scattered amplitude is then: 

A(q) = A 1 (q) + A 2 (q). 

A 2 (q) can be written in terms of Ai(q) as: 

A 2 (q) = A 1 (q)e incia e iqa/2 =A 1 (q)e iqa/2 , 

where a is the lattice parameter of the perfect lattice. The expected amplitude at Bragg 
position q=27i/a is then: 

A(27i/a) = Ai(27i/a) + A^ln/a)^ 71 = Ai(2n/a) - Ai(2n/a) = 0. 

For q values taken around 2n/a, the intensity is not 0, but follows the distribution 
displayed on the red curve of Supplementary Figure l.c. Two secondary maxima are 
found on both sides of q=27i/a, with an intensity that is half of the Bragg intensity 
associated to the perfect lattice. Moreover, the width of the two red peaks is the same as 
that of the Bragg reflection associated to the perfect lattice. These features are 
characteristic to the presence of a 71 phase defect. 

Dislocations are also n phase defects, but that develop in a three-dimensional 
lattice. Moreover, the phase variation can be smoother than the phase shift described in 
the case of the one-dimensional model. If one considers a perfect lattice 

R(x,y,z) = cos(27tx/a) + cos(27ty/a) + cos(27tz/a), 

the associated Bragg reflection will display the same features as the one presented for a 
perfect one-dimension lattice. Let us introduce a dislocation in that lattice, for example 
an edge dislocation with a line along z, and a Burgers vector along x. The distorted 
lattice can then be described as: 
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R(x,y,z) = cos(27ix/a+cp(x,y)) + cos(27ty/a) + cos(27iz/a), 



where cp(x,y) is the phase variation function accounting for the deformations induced by 
the dislocation. In the case of a dislocation line along z, with a Burgers vector along x, 
that function can be written as: 

cp(x,y) = -tt/2 * sgn(y-yo) + tan- 1 ((K/K x ) 1/2 *(x-x )/(y-yo)), 

where (xo,yo) are the coordinates of the dislocation line position, and K x and K y are the 
lattice force constants of the system along x and y respectively. This formula describes 
well the displacement field around the dislocation line, but not the atomic positions in 
the dislocation core itself. The lattice containing such a dislocation in the middle of the 
volume, with K x = K y = 1, and the associated phase variations are represented 
Supplementary Figure l.d and e. In that case, the phase variation is invariant around the 
dislocation line. Any cut of the phase distribution going through the dislocation line 
gives a similar phase profile as the one shown in red on Supplementary Figure l.b. The 
calculated Bragg reflection thus reveals the same features as in the one-dimensional 
case, namely a zero intensity value at Bragg position, and a torus of maxima around it, 
showing that the destructive interference takes place in every direction of space. The 
torus has an intensity that is half the intensity of the Bragg reflection associated to the 
perfect lattice, and the same width. The torus appears perpendicular to the dislocation 
line direction, and any cut of this intensity distribution by a 2D detector will reveal a 
split reflection, unless the cut is exactly perpendicular to the z* direction. 

The particular case of dissociated dislocation lines is discussed in the text. The 
principle is exactly the same, as the dissociated lines behave like a non-dissociated line 
in terms of phase around the lines. Only the stacking fault in-between the dissociated 
lines produces an additional signal. 
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Supplementary Figure 1 | Coherent diffraction pattern obtained when a n phase 
shift is present in real space, a, One-dimensional periodic lattice with no phase defect 
(blue line), and with a n phase shift in the middle (red line), b, phase associated to the 
lattices shown on a. c, Bragg reflection associated to the perfect (blue line) and phase 
shifted (red line) lattices shown on a. When a n phase shift is present, a destructive 
interference leads to a intensity equal to at the expected Bragg reflection position. The 
two maxima have same intensity, which is half the intensity of the Bragg reflection 
associated to the perfect lattice, and also the same width as it. d, Representation of a 
dislocation developping in a 2 dimensional elastically isotropic media, and e, associated 
phase, f, Intensity distribution around a Bragg reflection expected when the dislocation 
shown on d is present. Again, the intensity at Bragg position is exactly due to the 
destructive interference, and a torus of intensity appears when the media is elastically 
isotropic. 
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Supplementary Discussion 2 
Details about the intensity distribution seen on the 2D detector when 
the sample is slightly detuned from Bragg angle 



In this experiment, the illumination of the defining 7 um x 7 um aperture inserted 8 cm 
before the sample produces diffracted beams, impinging the sample with a slightly 
different angle than the direct beam, although the sample is still in the near-field region. 
The calculated beam received by the sample is represented on Supplementary Figure 2 
in linear and logarithmic scales perpendicularly to one set of blades. 
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Supplementary Figure 2 | Calculated beam profile at sample position. At a distance 
D = 8 cm from the 7 um x 7 um slits, the sample is positioned at the limit between the 
near-field and far-field regions. Due to the diffraction by the slits, the beam profile 
perpendicular to the slits blades already displays some fringes, which are supplementary 
incident beams impinging the sample, and are able to be diffracted if they enter the 
Darwin width of the reflection, a, Linear scale, b, Logarithmic scale. Fringes up to 100 
um from the main beam are still visible. 



The Ewald construction is represented on Supplementary Figure 3 for an incident angle 
equal to the (220) Bragg angle 6220- Supplementary Figure 3 shows that only the first 
order fringes appear in the vertical direction, but strongly attenuated. Indeed, the first 
order diffracted beam impinges the sample with an angle O220 +1.7 mdeg, which is still 
close to the (220) Bragg reflection condition (given the (111) monochromatized incident 
beam). The second order diffracted beam has an incident angle O220 + 3.5 mdeg, which 
no longer satisfies the Bragg condition. That is why it does not appear in that direction 
on the CCD detector. On the contrary, the fringes that are perpendicular to the 
diffraction plane are diffracted up to a high order, due to the reduced angular selectivity 
of the sample in that direction. 
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Supplementary Figure 3 | Image seen on the 2D detector when the sample is in the 
Bragg geometry. The sample is oriented in order to have the (220) reciprocal space 
node on the Ewald sphere. The supplementary beams due to the diffraction by the 
square slits defining the beam just before the sample are represented. These beams can 
be diffracted if they enter the Darwin width of the (220) reflection. The crystal 
truncation rod (CTR) decorates each spot in reciprocal space. The sketch only describes 
the behaviour in the diffraction plane, but the fringes and the associated CTRs are also 
present perpendicularly to the diffraction plane. Each spot or CTR that cut the Ewald 
sphere produces a signal on the 2D detector. A typical image obtained on the (220) 
Bragg reflection is shown. Only low order fringes appear in the diffraction plane, due 
to the small Darwin width of the (220) reflection. On the contrary, fringes appear up to 
high orders perpendicular to the diffraction plane, due to the lower selectivity in that 
direction. 
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Supplementary Figure 4 | Image seen on the 2D detector when the sample is 
slightly detuned from Bragg angle. The (220) reciprocal space node is no longer 
cutting the Ewald sphere. It only appears through the associated CTR that produces a 
high signal on the 2D detector (1). Some diffracted beams coming from the diffraction 
by the slits satisfy the Bragg condition and produce fringes with a higher counting level 
(2). In between, the CTRs associated to the other diffracted beams coming from the 
diffraction by the slits appear periodically between (1) and (2). Perpendicular to the 
diffraction plane, the CTRs associated to the other diffracted beams coming from the 
diffraction by the slits also appear. 
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